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ABSTRACT 


The purpose of this effort was to gain an insight into 
the Bi-phase cycle concept and to develop an experimental 
facility to test the key component of this power cycle--the 
Bi-phase nozzle. A system was developed and checked out that 
allows testing different Bi-phase nozzle shapes. The experi- 
mental system is provided with sufficient instrumentation to 
allow for the determination of the nozzle performance. In 
addition to the standard gas and liquid flow rates, measure- 
ment provisions exist for survey of the nozzle exit plane. 
Measurements can be made for the gas/liquid ratio and 


welocities at the nozzle exit. 
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I. INTRODUCTION 


Low speed prime mover systems in marine applications 
are attractive in that they negate requirements for heavy, 
voluminous and noisy reduction gearing. Direct drive con- 
ventional steam turbines and conventional electric drives 
have been tested and utilized with only limited success. 
What is required is a machine with low specific speed that 
is compatible at the same time with existing heat source 
systems. One option that may meet these requirements is 
the Bi-phase prime mover. The Bi-phase turbine is a 
machine having low specific speed and high torque. As 
such, it 1s ideally suited for marine propulsion. The Bi- 
phase turbine can drive the propellers directly without 
gear reduction. This feature leads to improvements in 
weight, volume and noise reduction. 

In addition, certain aspects of the Bi-phase concept 
may allow increased power production efficiencies. It may 
be possible to achieve near constant temperature expansion 
in the two-phase nozzle thus realizing a degree of 
"Carnotization" of the basic Rankine cycle. A detailed 
description of the Bi-phase power plant is given in Section III. 

A key element in the development of the Bi-phase prime 
mover is the Bi-phase nozzle. A Bi-phase nozzle is a flow 


accelerator in which fluid particles or droplets are 





accelerated by a vapor. The liquid and vapor may be 
chemically identical or dissimilar. The expanding vapor 
imparts kinetic energy to the liquid droplets. During 

this expansion, the droplets may also transfer heat to or 
from the vapor, resulting in a near isothermal expansion. 

In the nozzle, most of the available thermal energy of the 
vapor is converted to kinetic energy of the liquid droplets. 
The nozzle is converging-diverging, a requirement resulting 
from the low sonic velocity in two-phase mixtures. Computer 
codes [Ref. 1] exist that yield nozzle exit velocities to 
reasonable accuracies provided fluid properties are known, 
liquid particle size can be controlled, and one dimension- 
ally in the flow field 1s assumed. 

Tangren, Dodge, and Siefert derived the equations for 
isentropic flow of an immiscible two-component mixture 
with equal temperatures of the liquid and the gas. The 
first extensive experiments on two-phase nozzles were 
performed by Elliott [Ref. 1]. 

The purpose of this effort is to develop an experimental 
facility that will allow evaluation of two-phase nozzle 
concepts and associated liquid-vapor interaction phenomenon. 
The ultimate goal is to provide experimental capability to 
validate analytic models of Bi-phase nozzles. To this end 
an experimental system was designed that includes a high 


degree of flexibility to allow for rapid flow passage 





geometry variation. The instrumentation installed has the 
capability to sample liquid content at different points in 


the nozzle exit plane. 
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IL. BACKGROUND 


Turbine type prime movers, whether based on the Rankine 
or Brayton cycles, have many advantages. Among these are: 
high torque at low speed, low vibration levels as compared 
to internal combustion energy converters, and in some cases 
very high power densities. Practical considerations result 
in prime mover power outputs which are at relatively high 
rotational speeds. 

Power demand devices such as propellers and pump jets 
have rotational speed requirements set by considerations 
considerably different than those of the prime mover. Un- 
fortunately the rotational speeds of prime mover and power 
user seldom match. AS a conSequence, a power conditioning 
device, a gear box or reduction gearing is invariably 
mandated. In certain installations the reduction gearing 
may weigh more than the prime mover. 

From a purely thermodynamic consideration one would 
want the maximum cycle temperature to be as high as possi- 
ble. However, where the working substance must be 
contained by physical structures such as combustors, 
boilers, nozzles, and turbine blades, the cycle temperature 
limit is set by material considerations. Furthermore, 
secondary aspects of a prime mover, i.e., safety, may 


further restrict the range of cycle parameters. A case in 


ls 





point is the Navy nuclear power plant that does not admit 
Superheat. The turbine inlet conditions are those of 
Saturated vapor. This constrains the maximum cycle effi- 
clency to comparatively low values. 

A relatively new prime mover concept that may have the 
capability to circumvent some of the above identified diffi- 
culties is the Bi-phase prime mover. The technology base 
for the Bi-phase concept is founded in earlier work at the 
Jet Propulsion Laboratory on liquid metal power generating 
Concepts. 

The baseline Bi-phase concept is still the Rankine 
cycle (Fig. 1). The differences are subtle but of extreme 
importance. The working fluid is vaporized in the heat 
addition component of the cycle. This may be a conventional 
boiler generating high temperature and pressure steam. 
Normally the working fluid would then enter a nozzle/turbine 
system and energy would be converted to a mechanical form. 
In the Bi-phase system, however, a second fluid is intro- 
duced at boiler temperature and pressure just prior to the 
point where the steam enters the nozzles. This second 
fluid may have a vaporization temperature higher than water 
and thus will remain in the liguid state throughout the 
cycle. In the nozzle proper the vapor constituent acceler- 
ates and in so doing drags the droplets of the liquid phase 
to higher velocity (Fig. 2). This energy exchange of enthalpy 
drop of the vapor to kinetic energy of the liquid is the 


major departure from the conventional steam power plant. 
ile? 
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Downstream of the nozzle the high kinetic energy liquid 
and the expanded vapor are separated. The vapor is condensed 
and returned to the boiler. The ligquid is "parked" ona 
rotating wheel and energy 1s extracted by an impulse or pitot 
type turbine. After interaction with the hydraulic turbine, 
the fluid is returned for injection at the Bi-phase nozzle 
mixing section. 

It is the Bi-phase aspect of the system that is at the 
base of some of the potential advantages of the concept, 
which are: 

1. The energy imparted to the liquid phase is ina 
kinetic form of a high density mass. Hence, the 
liquid velocities involved are relatively low as 
compared to velocities of the vapor. Thus the 
output of the impulse hydraulic turbine will be at 
high torque and low rotational speeds--more in 
keeping with rotational speeds of driver elements. 

2. The energy extraction from the steam and the energy 
addition to the liquid takes place in an inherently 
small device--the Bi-phase nozzle. Thus, the high 
temperature, pressure, and velocities are confined 
to a high energy density component. 

3. Perhaps the most important aspect is the phenomenon 
that takes place in the nozzle proper between the 
liguid and vapor and how this interaction impacts 


the overall cycle performance. While it is normal 
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practice to assume "isentropic" expansion ina 
nozZle, it is not likely that this is what happens 
in a Bi-phase expansion. In all probability thermal 
energy is exchanged between the liquid and vapor 
constituents. This, in fact, may be the most 
attractive aspect of the Bi-phase concept. If 
thermal energy can be added to the vapor phase 

from the liquid, then the expansion can be made to 
approach an "isothermal" process. The net result 
is a degree of "Carnotization" of the basic Rankine 
cycle that can only lead to increased overall cycle 


efficiencies. 


Application of two-phase turbines to Navy propulsion has 


been examined in two design studies, References 2 and 3. 


In both cases, the following potential advantages were found 


when the Bi-phase turbine was compared to advanced steam 


turbine systems. 


jks 


Degiecreretency—-—Full load Output power performance 
Gams ranging from 20 to 50 percent were found. 
Direct drive--Direct drive at speeds ranging from 
90-4500 rpm was found possible with a single-stage 
turbine. 

Reduced volume--Volume reductions of 30 percent 
were estimated. 

Pre der lOg@eerrtcLlency—--Variable mass ratio 
enabled part-load (cruise) efficiency gains of as 
much as 100 percent. 
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An independent Navy study of a large Bi-phase marine 
turbine arrived at these same conclusions. These potential 
advantages lead to many possible applications for marine 
use. Most prominent among these advantages are those 
applications where space and weight are the most important 
considerations, for example, torpedoes, submarines, SES, 
SWATH, etc. 

It is thus apparent that the key element of a Bi-phase 
power system is the Bi-phase nozzle. It is this component 
that is the subject of this thesis. Specifically, the goal 
is to develop a research capability to study the flowfield 


and energy interaction in a Bi-phase nozzle. 
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Tt. Seer taoe CYCLE THERMODYNAMICS 


The Bi-phase cycle is basically a variation of the con- 
ventional Rankine steam plant. In reality, the "Bi-phase 
aspect" should be a multicomponent cycle as there is no 
requirement that the components (vapor and liquid) of the 
working substance be chemically identical. 

In a conventional Rankine cycle (Fig. 1) the liquid 
phase is introduced to the heat additive component of the 
system at the maximum cycle pressure. Energy is added and 
the substance is vaporized. This is a constant pressure 
energy addition process and follows a constant pressure line 
on the T-S coordinate system. The maximum cycle temperature 
may be well into the superheat region. Expansion, often 
idealized as isentropic, follows. During this pressure drop 
from boiler to condenser pressures, work is extracted by a 
turbine via a volume expansion. Constant pressure heat 
rejection returns the working substance to a liquid phase 
which then is pumped up to boiler pressure and reintroduced 
to the heat addition component. 

A common way of depicting the maximum cycle efficiency 
for any cycle is to apply the Carnot efficiency criteria. 
This Carnot efficiency depends only on temperatures of heat 


addition and heat rejection. 
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ie am 
n Tp 
Where = Cannot eLLicilency 
T, = Condenser temperature (OK) 
T = Boiler temperature (aK) 


Thus, a cycle operating between source-sink temperatures of 
232°c and 64°C respectively has a maximum energy conversion 
of 33%. The efficiency for a Rankine cycle with the same 
temperature limits and saturated vapor turbine inlet condi- 
tions yields an ideal efficiency of 29%. This efficiency 
is based strictly on thermodynamic considerations and does 
not consider irreversibilities such as friction losses, 
temperature Bobs in heat exchangers, and hydrodynamic losses 
in the cycle mechanical components. Furthermore, from 
practical considerations, the last stages of the expander 
see a very wet working fluid = 75% quality. This results 
in erosion and impact losses leading to further decrease 

in the power generator effectiveness. 

It is the Bi-phase cycle variation that has the potential 
to circumvent the above identified difficulties of the 
Rankine plant. 

The Bi-phase cycle adhieves improved efficiencies because 
the process involved allows the development of a near Carnot 
cycle. Thus, the Bi-phase concept achieves a degree of 
Carnotization of the Rankine cycle. This Carnotization 


takes place in the Bi-phase nozzle process. While in a 
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Rankine plant the expansion process is idealized as isen- 
tropic; in the Bi-phase system, it can be caused to approach 
an isothermal process. 

There are basically two Bi-phase techniques whereby a 


Rankine cycle can be caused to approach the Carnot cycle. 


A. TWO COMPONENT SYSTEM 

One method to Carnotize a steam cycle is to add a second, 
chemically different component that remains a liquid through- 
out the cycle process. This component provides a source of 
heat for the expanding vapor producing a near-isothermal 
expansion. In Fig. 1 it can be seen that the temperature at 
the end of the. nozzle expansion is 218°C compared to 64°C 
which would occur in the Rankine cycle. This results in 
a cycle efficiency of 31% versus 29% for the Rankine cycle. 
This cycle was invented by D. G. Elliott of Cal Tech for 
liquid metal MHD. [Ref. 1] 

The two-component cycle uses a low vapor pressure liquid 
and a high vapor pressure liquid having different chemical 
makeup. Some fluid combinations which have been considered 
are steam-Krytox, steam-Caloria, steam-lead bismuth eutectic, 
and Dow-therm-Therminol. The basic advantages of a two- 
component engine are low rpm, high efficiency with no steam 
extraction and a more compact turbine. Disadvantages 
include larger heat exchangers and the complexity of handling 
two fluids. Figs. 3 and 4 depict two possible methods of 


utilizing the two-component approach. 
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D. SINGLE COMPONENT SYSTEM 

In this cycle the expansion process is ideally paralled 
to the saturated liguid line that outlines the left hand 
side of the enthalpy dome. In order to achieve this type 
of Carnotizing, heat has to be removed during expansion. 
While it may be possible to do that through the walls of 
the nozzle, a more conventional way is to extract steam in 
steps. That means that a scheme of pressure staging is 
required, where after each expansion stage steam is ex- 
tracted and used to preheat feed water by condensing all of 
the extracted steam. 

Mae single component system has the advantage of elimin- 
ating o1l or heat transfer fluid as the second component. 
The two-phase mixture in this case consists of steam and 
water. Typical volume ratios of water to steam are 10:1 
or less so this cycle produces somewhat higher values of 
turbine speed. 

In order to obtain high efficiency several expansion 
stages are used. The five stage system resulting in the 
temperature-entropy variation of Fig. 5 is shown schematic- 
ally in Fig. 6. Steam from the heat exchanger flows to the 
first stage nozzle. A portion of the steam is also used 
to heat the water flowing to the first stage nozzle. The 
mixture of primary steam and water 1S expanded in the 
nozzle from an inlet vapor quality of 15% to an exit 


quality of 18%. The two-phase mixture drives a rotary 
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Separator turbine. Part of the separated steam is extracted 
and flows to a feed water heater. The remainder of the 
steam is remixed with the water exiting the first stage 
turbine. This mixture is expanded in the second stage 
nozzle from an inlet vapor quality of 12% to an exit 

Suelity of 18%. 

The two-phase mixture drives the second stage rotary 
separator turbine, where additional steam extraction occurs. 
In the final stage, steam exhausting from the turbine flows 
Pememeecondenser. The condensate 1s remixed with the 
liquid exiting the last stage turbine. The feed water 
flows through several stages of a feed water heater where 
the extraction steam 1s condensed on the feed water. 

Five stages are shown in Fig. 6. The actual number of 
stages will be determined by a tradeoff between efficiency 


and cost and complexity. 
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TV. BisPaASE NOZZLE THEORY 


The Bi-phase nozzle is the key element in the Bi-phase 
engine concept. While initially it may appear that the 
design of the nozzle profile is quite straight forward, 
closer examination will reveal that the question is con- 
Siderably more complex. The complexity steams from the 
nature of the working fluid itself. It is a mixture of 
liquid particles and a compressible gas or vapor. 

The complicating factors can be grouped into three 
categories: 

1. The working fluid has two constituents significantly 
different in density. Thus, there is an aerodynamic 
interaction between the components (liquid droplets 
are accelerated by the vapor). 

2. The shape and size of the liquid particles is 
governed by the forces on them. In this particular 
case the surface tension and inertia forces are of 
key interest. 

3. Thermal interaction between the working fluid 
components can range from none at all to phase 
change of the liquid to its vapor form. 

Another factor of interest in Bi-phase flow is the 

question of the velocity of propagation of a disturbance, 


i.e., velocity of sound. While the velocity in a particular 
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vapor and in a particular liquid are well established, the 
Situation in a mixture of vapor and liquid is considerably 
different. 

Consider the following discussion from Olson [Ref. 4]. 
The acoustic velocity in any medium is c = /K/p, where K is 
the modulus of elasticity for the medium and o is the 
density. For gas-liquid mixtures the acoustic velocity 
becomes less than that for either the liquid or gas alone. 
For aliquid with a small concentration of gas nuclei, the 
elastic modulus of the mixture is reduced, with no appre- 
ciable reduction in density, and thus the acoustic velocity 
is reduced. For a gas with minute liquid droplets, the 
density of Pree oceans 1S increased with no appreciable 
change in elastic modulus, and again the acoustic velocity 
for the mixture is reduced. The velocity of sound in a 
14 percent water and 86 percent air mixture is about 30 m/s 
(depending on the pressure, temperature, and impressed 
frequency) as compared with about 1460 m/s for water and 
about 335 m/s for air. 

The flow phenomenon of a Bi-phase mixture has been 
analyzed in Ref. 1. It is repeated as follows. The 
problem is illustrated in Fig. 7. A spatially uniform two- 
component mixture of liquid drops and gas enters a nozzle 
at high pressure and low velocity and expands to low 
pressure and high velocity. The objective of the analysis 


is to determine, for a specified pressure, the drop diameter 
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L’ velocities Ve and Vy, and 


flow rates mM, and m. of the gas and liquid phases, respec- 


D and the temperatures be and T 
L 
tively, at each station inthe nozzle given the initial 


values of D, T IL V V the total flow rate, and the 


on ig re ag 
properties of the fluids. 
The five relations employed to compute the five unknowns 


and V are (1) the momentum equation for the 


i eg 
mixture, (2) the energy equation for the mixture, (3) the 
drop drag equation, (4) the drop heat transfer equation, and 
(5) the drop breakup criterion. Solubility and vapor 


pressure relations provide the flow rate ratio m,/ My, « 


A. ASSUMPTIONS 

The assumptions employed in the two-component analysis 

are as follows: 

1. The liguid is uniformly dispersed as spherical drops 
all of the same diameter. 

2. The drops break up to limit the Weber number to 6. 

3. There are no external forces acting on the two-phase 
mixture other than pressure and wall shear, and 
there is no heat transfer to or from the mixture. 

4. The flow is one-dimensional. 

5. The drops are large enough for the surface curvature 
to have negligible effect on the vapor pressure of 
the liquid and for the surface energy to be negligible. 

6. The drops are isothermal. 


7. The gas mixture obeys the additive-pressure law. 
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8. The partial pressure of the predominantly liquid 

component in the liguid is given by Raoult's Law. 

9. The concentration of the predominantly gaseous 

component in the liquid is given by Henry's Law. 
10. The volume of the liquid solution is equal to the 
sum of the volumes of the pure liquids. 

Assumption 1 restricts the analysis to nozzles having 
spatially uniform injection of the liquid into the gas and 
operating at gas-to-liquid volume ratios greater than 
unity. Assumption 2, the drop breakup criterion, states 


that drop diameter is limited to a value D for which ile, 


2 
Ge = MeeD/2Z20 = 6. Thus 
B Pag 
a eo 
vee . 2 (1) 
Oo V 
9g 


where Og is the gas density, V5 is the slip velocity ep a 
and o is the liquid surface tension. The form of Eq. (1) is 
physically reasonable in that the Weber number is proportional 
to the ratio of stagnation pressure 0 Ws /2 to surface 

tension pressure 40/D. Hence, a drop would be expected to 
flatten and break up at a sufficiently high value of We. 

This has been verified experimentally and the critical Weber 
number found to be 6, within a factor of about two. An 
additional restriction is that for actual breakup to occur, 


the time spent at a Weber number exceeding 6 must be longer 


than the natural period of oscillation of the drop, 
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1 (p,D9/) 2/4, where Oy is the density of the liquid. This 
requirement is met only in two-phase nozzles longer than 
about 10 in. and Assumption 2 may cause the analysis to 
overestimate the exit velocity by increasing amounts as 
the nozzle length decreases below 10 in. 

Assumption 3 excludes magnetohydrodynamic and mechanical 
body forces. The exclusion of wall heat transfer 1s correct 
for the insulated nozzles of interest for power systems. 

In addition, the relatively high velocity results in short 
residence lines in the nozzle proper. 

Assumption 4 is closely met in practical nozzles since 
good performance requires small wall angles, large throat 
radius of curvature, and uniformly distributed injection 
of the fluids at the nozzle entrance. 

Assumption 5 is valid for the drop sizes of 0.001 to 
0.010 in. produced by the Eq. (1) breakup criterion. 
Assumption 6 is valid because of the rapid internal circula- 
tion in drops. Assumption 7 introduces negligible error 
in most cases of practical interest since the vapor pressure 
of the liquid is small and needs only to be evaluated 
approximately. 

Assumptions 8, 9, and 10 are either valid, or cause 
little error, for fluids of low miscibility, which are the 


fluids of interest. 
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B. DERIVATION OF EQUATIONS FOR FREE-STREAM FLOW 
i) CeOnt inul ay 
Referring to Fig. 1, the nozzle flow area A is equal 


to the gas flow area m,/0 V_ plus the liquid flow area 








Gg 
m,/P,Vz- Paws, 
IL Ds 
A = mM + Zz) 
g Pg PV) 


where r is the mass mixture ratio my, /M,- 
2. Momentum 
By Assumption 3, the only force acting on the free- 
Stream flow is that due to the pressure gradient. If M is 
the momentum flux at flow area A, the change in momentum 


flux across pressure increment dp is 
aM = -Adp (3) 


The momentum flux can be written as the sum of the 
momentum fluxes of the gas and liquid. MThus, 


M = mVg + mV} (4) 


If the flow were allowed to continue at constant 
pressure, Vy and Vy, would become equal to each other at the 
mass-weighted mean velocity V. Since for this process, 

dM = 0, the value of V is given by 


(m., a m,)V = ee a5 mV, (5) 
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or 
- Vests rv, 
ae (6) 


Thus, the momentum flux can be written 
M = a) + m,)V (7) 
Since m, 2 m, 1s constant, the change in momentum flux is 


dM = (m + m.)daV (8) 


Substituting Eqs. (8) and (2) into Eq. (3), dV can be written 





= aly ul Me 
av = ———/{ —— + dp (9) 
l+er (sa, =) 


The slip ratio is defined as 


V0 
ll 
< 
—< 
< 

| 


Ng - V,)/V ‘Gly 


This equation can be combined with Eq. (6) to give Ve and 


V. in terms of V: 


L 
‘e = ( + 72.) Vy any (leis) 





V bV (Cl) 


or 


| 
> aie. 
1 
a 
+10 
ry 
ae | 
<i 
il 


The gas density can be expressed as 


= W RT 13 
ae ae 5 (eleS) 


where We is the effective molecular weight of the gas mixture 
and R is the universal gas constant. Eq. (13) 1s the 


definition of the effective molecular Megs Nig which is the 
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quantity that gives the actual gas density when substituted 
moma. (13). 

Substituting Eqs. (11)-(13) into Eq. (9), the 
differential momentum equation is 


RT 
SS =P. 2 g + © (14) 


The quantities a and b are slowly varying because s is 
typically only 0.1 to 0.3 and slowly varying. The quanti- 
ties r, Tee i and op, are also slowly varying. Integrating 
Eq. (14) over a pressure increment Ap, for which a, b, r, 
a W_, and Pp, are constant to within the desired accuracy, 


g 
the change in v7 is 


Ap 
= lee a 
AV = - iy awWP + bo, dp (15) 
o- ap 
2 


All quantities other than pressure can be taken outside the 
integral and evaluated at their mean values (denoted by 


subscript m) corresponding to the mid-interval pressure p. 





Thus, 
- 2 
AV2 = - 
1+ rn 
A A 
p+ = pt = 
RT q Cn 
x 4m a Ee dp }(16) 
m m*” 1 
Sin = Ap m p- 4p 
2 2 
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Performing the integrations, 





RT 
g , 1g EA, 
we > = ae sur lise, & : 7 55 amie |) (LD 
m mg P P m°~m 


Equation (17) is the final form of the momentum equation. 


oe Snergy 


The enthalpy change of the mixture between state l 
(the beginning of pressure interval Ap) and state 2 (the 
end of the interval) can be evaluated in two steps: 


(1) phase change at Pye AB T and (2) change to Po, T 


eye ical L. 


ie at £Lixed composition. 


g2 
The enthalpy change for step 1 is 


enthalpy required 


, 
AH. = amount of A, x to vaporize and 
HL vaporized | heat unit mass of 
A from Ty Ee) Tot 
enthalpy required 
amount of B x to vaporize and 
vaporized heat unit mass of 
B from Vr, Gort 
1 J 
kinetic energy 
— oe required to 
S =e x accelerate unit 
P mass from Vr to 
V it 
Fy 
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or 


AH, = (m Sei) oe etc Ze = PT. )] 
| 3. 7; ay oh Jy 1, 
+ (m = Vika. WES (T -T. )] 
a “By, by Ba S) L, 
+ (h -m_)(vV% = v2 )/2 (18) 


J5 i hl 1, 


where L and c are latent heat and specific heat, respectively. 


Introducing more compact notation, 


_ JG + 
+ he + Cy 6,7) + 5 (19) 


The enthalpy change for step 2 is evaluated from the 
temperature, pressure, and velocity changes, with properties 


evaluated at mean T and p for the interval. 


_ 1 2 2 
P= seme or 4) + © ee - v2 4) 
2-90 Gm GQ CQ A GQ 
Eo PY 
+m, fc (oe - T. ) + 
1, tn fs 1, Py 
m 
i 5 (v2 - v2 )] (20) 
2 i 


oy, 





= m Ge Az ee) 
Go\ 9, 9 2 
2 
AV 
fu my Cy AT, + Ap + -=* (21) 
2 m el, 


By Assumption 3, no work is done by the free-stream flow and 


no heat is transferred to it. Hence, 


AH, ct AH, = 0 (22) 


Suestitucing Eqs. (19) and (21) into Eq. (22) and 


solving for AT_ gives the energy equation for the mixture: 


i Ave r ave 
A eo eo 
g e. 2 2 i 1 

Im a pl 
ams, Mag 

+ ah + a + Cc 6,7) 

22 32 oii 

Am, 

+ = ity, + Cy 50) (23) 

<5) an 
4. Drag 


Although no force other than pressure acts on the 
free-stream flow as a whole, a drag force exists between 
the phases. Hence, a second momentum equation must be 
written uSing as the control volume the boundary between 
the phases. 

The two forces acting on each liquid drop are the 


buoyancy due to the pressure gradient and the drag due to 


38 





the relative gas velocity. The sum of these is equal to 
the mass times the acceleration of the drop. Thus, for a 


Single drop. 


dynamic pressure of ~ (s hat=\o Mole has frontal area 
relative gas flow efficient of drop 
_ | volume XY pressure _ mass accelera- 

of drop gradient of drop Elon Onmadtonm 


or 


« og l¥ V2) Cy mD* | mD2 dp 
4 


= mo, \(y oe (24) 
6 "L/\*1 dx, 


The absolute value sign in the first term makes the 


drag force positive when Ye ee and negative when west 


L 
poerevinG EG. (24) Tor dv 


L° 


Tin? 
So_|s|sV2c_dx 
BS O- _ _ dp 
avs _ AO VvAD V (25) 
alee aa 


DPerrerentiating Eq. (12), dv, can also be expressed in terms 


Sees, c;, ana V. Thus, 





_ = sdr _ ds } (26) 
Simo” * Vale Se Tapes Z 


Solving f£&or ds, 


a ( 
gg = PULtridv, _sdr_ (27) 


V 


| 





Suse recut tig dv, Baemleng. (25), noting that dv = ave /27, 
using Eq. (12), and writing for a finite increment, results in, 
b (1 +r_)AV2 (1+ r_)Ap s_Ar 
ss = SH Sod 
2Ve DAP L Mies m 
m 
30g [Sql SQCp_ (4 + rx) Ax 
- EE (28) 


4b P41 D 
m 
This is the drag equation employed when x is specified 
mea function of p. 
Solving Eq. (28) for Ax yields the required alterna- 


tive equation: 


b2o, Av2 
AD a Ln 
Ax = 5 RYO) Gs 
V 2 
30g [Sn [SnCp V 
<2 
b Pym / Smor ] 
= | 
Vices mer 9/1 ‘+7? 
m m 


5. Heat Transfer 
Although no heat is transferred to the mixture as 
a whole, heat transfer exists between the phases. Hence, 
a second energy equation must be written using as the 
control volume the boundary between the phases. 
The work dW done on the liquid is that due to drag 
by the gas. (Only work done by shear or shaft forces is 


included in @W when writing the First Law for a control 
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volume). Multiplying Eq. (24) by the number flow rate of 


drops N = eit, /mD°p,, the drag force F. on that quantity of 


d 


liquid is 


m av 
= 2. 1 op i 
Fa gg IVgIVEpmD , ax + mV, a (30) 
The work done on the liguid is 
av? 
-dW = F.dx =m (dp + | (ot) 


The heat dQ transferred from the liquid is made up 
of two parts: (1) the convective cooling due to the 
temperature difference between the liguid and gas and (2) 
the evaporative cooling due to the latent heat supplied to 


the liguid vaporized. The convective cooling is 


sdOee> ha N(T, - T,) dt (32) 


d 


where h is the heat-transfer coefficient, As = nD? is the 


surface area of a drop, and dt = dx/V_ is the time required 
to traverse dx. Thus, 
6hm, (T_ - T,)dx 
-aZ = ae) Ok: en (a) 


C Do 
7a 


The evaporative cooling is 


-dQ., = Hama + ee (34) 
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The change in enthalpy of the liquid over the pressure 


increment dp is 


° ( dp avy 
dH = mM c, aT, + = + ——— (35) 


\ i_ 
Subctitucing Hos. (51), (33), (34), and (35) into the 


steady-flow energy equation dQ - dW = dH, the result is 


6hm, 6Tdx 


Dev, = Ldm, = “poy | = m.c,daT (36) 


Hie lee 1. 
S| 


where 6T = T -T._. 
g L 


Writing for a finite interval, the final form of the 


drop heat-transfer equation is 


e p a 
; éhé_TAx vai Am, 
oo = pa sae = | ee (37) 
IL fe: Do, V a mf bom 
Ie 1, 1. m 1.4 m aoa 


Equations (1), (17), (23), (28) and (37) are the five 
equations that must be solved simultaneously to obtain the 


V and 


Te” g’ 
V, as a function of the independent variable p. To carry 


out the solution all quantities in the equations must be 


values of the five dependent variables D, Ua z 


expressed in terms of these six variables. 

The preceding equations form the basis for the 
mathematical model which is used to predict, based on 
inlet conditions, the exit velocity, and temperature of the 
mixture. These equations also form the basis for the model 


which provides the optimum nozzle shape given a set of inlet 
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conditions. Some additional relationships are, however, 
required. These are: 

1. Phase properties--to establish the mass ratio, 
mass flow rate ratio of gas to liquid, and the 
thermal conductivity of the mixture. 

2. Liquid drop drag coefficients. 

3. The liquid drop heat transfer coefficients. 

4. Boundary layer momentum thickness and displace- 
ment thickness. 

Ta omen Lrictilon coefficient. 

These five additional relationships are developed in detail 
i set. 1. 

The preceding "basic" model in conjunction with the 
additional equations form the core for four computer programs 
which provide real and isentropic flow output parameters, 
and an optimized nozzle contour given a set of inlet 
conditions. 

Reference 1 reported that the optimum nozzle profile 
has an elongated throat region, that the required shape is 
insensitive to drop diameter and nozzle length, and that the 
exit velocity is insensitive to departures from the optimum 
shape. For given fluids and pressures, the exit velocity 
was found to be insensitive to nozzle length and flow rate 
beyond certain minimum values. This insensitivity of 
nozzle exit velocity to nozzle shape and size results from 


the inverse dependence of drop size on slip velocity, which 
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apparently acts to adjust the slip to a value that holds the 
exit velocity in the range 75 to 90% of isentropic. Thus, 
it is difficult to design either a very good or a very poor 


two-phase nozzle. 
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V. EXPERIMENTAL APPARATUS 


The experimental apparatus was deSigned to provide a 
facility to study the performance of a Bi-phase nozzle. 

The basic requirements were to allow overall nozzle per- 
formance evaluation as well as optical diagnostic probing 
of the Bi-phase flow-field. 

The experimental system can be conveniently grouped into 
Me~esetosyscems. These are: the nozzle proper, nozzle 
assembly, air supply system, liquid injection system, and 
nozzle exhaust plane instrumentation. Each subsystem is 
described in the following sections. Figure 8 is an overall 


system drawing. 


A. NOZZLE 

The nozzle was designed based on recommendations given 
in Refs. 1 and 5. It is convergent-divergent, 12 in. long 
with a throat area of 0.25 ine Inlet area is 2.0 anaes 
The throat is located 4 in. from theinlet (1/3 of the over- 
all length). The nozzle is constructed by "sandwiching" 
two & inch thick machined aluminum nozzle profile plates 
between }; inch plexiglass plates (Figs. 9 and 10). The 
aluminum nozzle plates are located at the end of a 30 inch 
long "test" section and are easily replaceable. The test 


section provides the necessary space for pressure taps, 


liquid injection tubes, and other instrumentation and 
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—==—s ie NOZZLE ASSEMBLY 


Inlet Flange 


Plexiglas 





Aluminum Plate 


Fig. 9. Nozzle Assembly Drawing. 
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Nozzle Assembly Photograph. 
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configuration options. The ability to replace the nozzle 
profile plates allows for modification of the nozzle 


configuration. 


B. AIR SYSTEM (Fig. 11) 


mom? £t.> 


tanks supplied with compressed air from an 
Ingersol Rand three cylinder air compressor provide the air 
storage volume required to support nozzle operation. The 
nozzle is supplied with air at variable pressure via 3" ID 
piping. The air supply piping contains two pressure control 
valves, a solenoid actuated nitrogen operated 3 in. ball 
valve, and a flow measurement orifice plate. The orifice 


flow coefficient was determined from Ref. 6, Table 4. The 


air mass flow rate can be determined from the equation 


W. = 359 K d“Fa Y Vay 


Wr - air mass flow rate in LBm/HR 


K - orifice flow coefficient Table 4 [Ref. 6] 
ad - orifice diameter 


F_- thermal expansion factor Fig. 38 [Ref. 6] 


- expansion factor, Fig. 40b [Ref. 6] 


- differential pressure in H.~0 at 68°F 


2 
Y - Specific weight of the flowing fluid at the inlet 


side of the primary element (LB/FT?) 
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and 


a Sr ee 
W 627 
where 
Bes - inches of differential of inches of manometer fluid 
Ym 7 specific weight of manometer fluid LB/FT> 
Y, ~ specific weight of fluid separating manometer from 
flowing fluid 
h- differential head, in. of H,0 at 68°F 
And 
y Se 
Where 


le Specific weight of dry air at actual temperature. 


For an example calculation, see Appendix 2, example 5 from 


Ref. 6. 


Geese olD INJECTION SYSTEM (Fig. 12) 

The liquid injection system consists of a pressure 
vessel to contain the water which is pressurized from a 
nitrogen tank. The pressure vessel supplies pressurized 
water to the injection tube via a solenoid stop valve, a 
flow control valve and a flowmeter. The flowmeter was 
calibrated using an accurate scale and a stop watch. 

The calibration data and chart are recorded in Fig. 13. 


The liquid injector is a 0.25 in. brass tube inserted in 
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the 3 in. ID air supply pipe just upstream of the flange 
connection to the test section. The injector tube is 
drilled with 5 holes facing the test section entrance. 

It has been proposed that more efficient liquid atomization 
could improve nozZle efficiency; therefore, the drilled 
holes were made as small as possible consistent with 


achieving a significant liquid mass flow rate. 


D. NOZZLE EXHAUST INSTRUMENTATION 

A cylindrical drag body, which is attached to a rectan- 
gular cantilever beam, is placed in the nozzle exhaust 
flow stream. The drag body unit is attached to a bracket 
which is attached to a rectangular transversing plate 
(Figs. 14 and 15). The transversing plate is screw driven 
to move in a transverse direction across the nozzle exhaust. 
This movement causes the cylinder to move through the 
entire nozzle exhaust flow stream. By placing strain gages 
on the drag body unit, the cantilever beam deflection can 
be determined and related to the force on the cylindrical 
portion of the drag body. That force, in turn, can be 
related to the local velocity of the flow stream and these 
velocities plotted to give a velocity distribution at the 
nozzle exhaust. A similar plot can be constructed for the 
momentum distribution 


Z 
qa ou 


Drag (LB -) = Ch 5 
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Instrumentation Photograph. 
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where 


C. - Drag coefficient of the cylinder 


- (Smooth cylinder 1/8 in. in diameter) 


A - Projected area of the cylinder 1/2" XK 1/8" = 1/16 in. 
0 - Density of the mixture 
U - Local velocity of the exhaust flow stream 


Selving for U: 


(iD aera 
Pee Dao LOrce 
le Ty pC nA 


Figures 16 and 17 are calibration data for the drag cylinder. 
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Vis - GONCLUSIONS 


Successful Bi-phase nozzle operation was achieved and 
Sustained under controlled conditions. Liquid injection, 
nozzle, and instrumentation system all operate satisfactori- 
ly. The liquid breakup and acceleration in the nozzle 
passage was observed. 

With the nozzle pressures employed, i.e., 20, 30, 40, 
and 50 psig an instability was observed in the downstream 
portion of the nozzle passage. It was apparent that, for 
the nozzle geometry installed, there was a transition from 
supersonic le in the throat region to subsonic conditions 
in the diverging downstream passage. It 1s apparent that 
at this point the passage converts from a supersonic nozzle 
to a subsonic diffuser. The transition point appeared to 
be unsteady, initiating alternatively from the sides of the 
nozzle. The position at which transition takes place could 
be controlled at will. As the liquid mass ratio and/or 
upstream pressure were increased, the passage tended to 
perform aS a supersonic nozzle, i.e., the supersonic- 
subsonic transition was driven further downstream towards 
the exit. This is as one would expect. As the liquid 
mass ratio iS increased, the sonic velocity decreases and 
the passage tends to a supersonic nozzle. Also, as the 


pressure ratio increases, the shock or transition to subsonic 
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conditions is driven downstream. An exit plane survey 

(Fig. 18) using the instrumented drag cylinder confirms 

the existance of a transition point somewhere in the diverg- 
ing portion of the flow passage. Thus, one is led to 
conclude that the role of change of flow passage area is 

too high for the test pressure ratios. From Fig. 18 it 

is also evident that the flow tended to oscillate from one 
side of the nozzle to the other. This was confirmed by 
visual observations as well as audible pulsations. 

It would be desirable to develop a degree of flexibility 
in changing or altering the nozzle profiles, particularly 
the diverging passage. Originally it was envisioned to 
replace the nozzle inserts; however, it may be more expedi- 
tious to provide an infinitely and continuously variable 
system. This may be achieved by a flexure system controlled 
by jack screws in the exit plane. The second desirable 
aspect would be the capability to diagnostically probe the 
flow field in the passage proper. A laser anemometer may 
provide the capability; however, no information exists on 


the use of LDA in Bi-phase flow. 
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